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METHOD AND APPARATUS FOR BLOOD VESSEL 
PARAMETER DETERMINATIONS 


BACKGROUND OF THE INVENTION 
FIELD OF IHE INVENTION 
The present invention relates to blood flow diagnosis in general, and to 
a method and ^stem for determining the coronaiy reserve and the relative 
coronary reserve, in particular. 

DISCUSSION OF THE RELATED ART 

The coronary velociiy reserve (CVR) also known as the coronary flow reserve 

(CFR), and generally referenced hereinafl;^ as coronaiy reserve, is defined to be 

the ratio between llie maximal and the resting coronaiy blood flow. The maximal 

coronaiy blood velocity occurs at times of highly intensive exercise, or it is 

induced when the subject is injected with coronary vasodilator such as adenosine. 

Tlie coronaiy reserve represents the ability of the coronary arteries to supply the 

excess blood needed to comply with the excessive pumpmg requirements, 

relatively to normal (resting) conditions. The normal coronaiy reserve is greater 

than 3.0 and in some individuals is greater than 5.0. When the CFR is impaired 

(values equal to or less than 2.5) the coironaiy arteries are unable to supply the 

excess required amounts of blood. Reduced CFR is often associated with angina 

pectoris, diabetes mellitus, systemic sclerosis. Coronary Syndrome X and other 

clinical conditions. Determining the CFR is also valuable when assessing the 

severity of a stenosis in a coronary artery, smce there are cases when a stenosis is 

found, but is not the sole cause of ischemia. Therefore, prior to performing a risky 

and expensive medical operation, determining the CFR should be considered. 

Flow measurement is also useful for predicting long-term success of treatment 

and comparison of efficacy of various treatments. Some currently available 

techniques for measuring the CFR include invasive methods such as a Doppler 
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catheter, or a pressure wire, which follows a leading catheter that emits tibie 

contrast agent The Doppler catheter or pressure wire methods have several 

disadvantages* Extra intervention procedures are required, which makes the 

Doppler catheter or the pressure wire methods more costly, more invasive, and 

therefore more risky. In addition, the extra wire or catheter affect the flow itself 

and therefore impair the measured flow. Yet another drawback of the Doppler 

catheter and the pressure wire is that the catheter or wire must be accurately 

aligned, otherwise the results are impaired. Another method for measuring the 

CFR involves Magnetic Resonance (MR). This procedure is used for diagnosis 

only while the catheterization procedures are used for therapeutic purposes as 

well- The MR method is also costiy and physicians are reluctant to use it 

Another method for measuring the CFR is the Digital Subtraction Analysis 

(DSA). When using the DSA method, the flow is evaluated directiy firom the 

angiogiams, by measuring the change in the gray level in a specified region, 

caused by the spreading of the contrast agent. This method uses the relation 

between the gray level and the volume of the contrast agent (as in densitometry) 

to calculate the flow directiy from the change in gray level between successive 

frames. A limitation of the DSA methods is the difficulty in calibrating the 

relation between the gray level and the volume. This relation depends on imaging 

conditions such as the X-ray energy, the extent of magnification, the distance 

between the source image and the receptor (SID), or others, and contrast material 

parameters. There are too many parameters and not enough available information 

to accurately determine the CFR. Another drawback of the DSA method 

originates firom the area captured in the x-ray. Since this area is substantially 

larger than the area of the artery, the error calculations are mush larger, which 

again harms the measur^ents. 

Yet another method for evaluating the CFR is the Contrast Propagation 

Algorithm (CPA). When using the CPA method, the flow is determined by 

observing the propagation of the haze. The limitations of the CPA are that it is a 

difiScult task to follow the haze, due to its indistinct nature, and that it is required 
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to know the exact structure of the artery in order to take into account the 3D 
geometty. Therefore, this method is hard to execute and suffers form inaccuracies. 

A parameter related to the coronaiy flow is the TIMI Flow Grading, which is 
a qualitative assessment of dye washout during contrast angiography. Usiag the 
TIMI flow grading for assessing the coronary reserve is qualitative and subjective, 
and therefore does not provide accurate measurements. 

There is therefore a need for a relatively cost effective, easy to use, accurate 
and reliable method and system for measuring the C3FR of a subject and other flow 
related measurements. It is desirable that the melhod will provide accurate results 
for flow measurement at different times during the cycle of the heart beat, and 
will incorporate, as a preferred embodiment, the full heart beat cycle for providing 
accurate average results. It is also desirable that the method can be used during the 
catheterization, or at a later tune, and will not imply extra invasiveness beyond a 
standard catheterization. 
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SUMMARY OF THE PRESENT I^^S/^^ 
It is an object of tiie present invention to provide a novel method and 
apparatus for the determination of coronary reserve and relative coronary reserve 
of a specific coronary artery segment and other flow related measurements. 
In accordance with the present invention, there is thus provided a method for 
determining the arterial reserve of a subject having a blood flow of a velocity 
below or equal to a maximal velocity value, the method comprising the following 
steps: receiving one or more first models, each model representing one or more 
substantially fioced segments of one or more arteries of the subject at a plurality of 
points in time associated with and during one or more parts of one or more first 
heart beat cycles, the artery segment having a proximal cross section and a distal 
cross section; injecting contrast agent into the artery of the subject, when the 
subject is in a non-hyperemic state, the contrast agent is injected at an injection 
area having a proximal cross section and having a distance firom the substantially 
fixed segment; determining one or more first parameters from one or more first 
angiograms representing the substantially fixed artery segment, said first 
angiogram taken form a projection angle; injecting the subject with substance that 
simulates hyperemia; receiving one or more second models, representing the 
substantially fixed segment of the artery of the subject at a plurality of points in 
time associated with and during one or more parts of one or mo^e second heart 
beat cycles; injecting the contrast agent to the artery of the subject, said subject 
being in a hyperemic state, said contrast agent is injected at the injection area; 
determining one or more second parameters fi-om one or more second angiograms 
representing the substantially fixed segment, said second angiogram taken form a 
projection angle; and determining the arterial reserve as the ratio between one of 
the first parameters or a combination thereof and one of the second parameters or 
a combination thereof. The first parameter can be a ratio between the maximal 
velocity value of blood within the artery segment when the subject is at a 
hyperemic state and the distance between the substantially fixed segment and the 
injection area. The second parameter can be a ratio between the maximal velocity 
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value of blood within the arteiy segment when the subject is at a non-hyperemic 
state and the distance of the aiteiy segment from the injection area. Within the 
method, deteraiination of the first parameters can comprise the steps of: 
detemiining from the first angiogram tsken at a predetermined projection angle 
and the first model, a density curve for the artery segment; obtaiiiing a velocity 
profile of the blood flow within tiie artery; performing curve fitting for the density 
curve to determine the first parameter. Within the method, determination of the 
second parameter can comprise the steps of: determining from the second 
angiogram taken at a predetermined projection angle and second model, a density 
curve for the at least one artery segment; obtaining a velocity profile of tfie blood 
flow within the artery; performing curve fitting for the density curve to determine 
the second parameter. The distance between the injection area and the 
substantially fixed segment can be the distance between the distal ctoss section of 
the injection area and a cross section of the sfubstantially fixed artery segment 
located , at equal distances from flie proximal cross section and from tihie distal 
cross section of the substantially fixed artery segment. The method can further 
comprise the step of creating the first or the second models of the artery of the 
subject. The first or the second models can be three-dimensional models. The 
method can fiirther comprise the step of determining the projection angle and the 
voliraies of the fixed segment of the arteiy of the subject at a plurality of points in 
time associated with and during the one or more parts of the first heart beat q^cle, 
from the first model and the stq> of determining, the projection angle and the 
volumes of the fixed segment of the at least one arteiy of the subject at a plurality 
of points in tinie associated with and during the one or more parts of the second 
heart beat cycle, from the second model. The method can ftirther comprise the 
step of compensating for the non-perpendicularity of the substantially fixed 
segment of the artery of the subject The method can fiirther comprise the stejp of 
registering the first angiogram with the first model or flie second angiogram with 
the second model. The method can fiirther comprise a step of determining TIMI 
grades from local gray level curves in multiple points of the arteiy. The method 
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can further comprise the step of detenzuning a relative arterial reserve as the ratio 

between the arterial reserve determined for a first arteiy segment and the arterial 

reserve determined for a second arteiy segment The first arteiy segment can be 

diseased or suspect as being diseased and the second arteiy segment can be 

healthy. The arterial reserve can be arterial coronary reserve. The contrast agent 

injection can performed during the systole of the subject, or continuously 

throughout an integer number of heart beat cycles of tiie subject. The contrast 

agent can be injected radially. 

Another aspect of the present invention relates to a method for determining die 

blood flow output of a subject having a blood flow having velocity values below 

or equal to a maximal velocity, the method comprising the steps of: receiving one 

or more models, representing one or more substantially fixed segments of one or 

more arteries of the subject at a plurality of points in time associated with and 

during one or more parts of one or more heart beat cycles; injecting contrast agent 

into the arteiy of the subject at an injection area having a distance firom the 

substantially fixed segment; determining firom one or more angiograms taken at a 

projection angle and the model, a density curve for the substantially fixed 

segment; obtaining a velocity profile of the blood flow within the substantially 

fixed segmeait; performing curve fitting for the density curve to determine one or 

more parameters; substituting the at parameters in the . velocity profile to 

determine velocity values; and integrating the velocity values over a cross section 

of die substantially fixed segment to obtain the arterial segment output values of 

the blood flow within the substantially fixed segment of the artery. Within the 

method, the parameter can be the ratio between the maximal velocity of blood 

within the substantially fixed segment and the distance of the substantially fixed 

segment from the injection point Within the method, the distance between the 

injection area and the substantially fixed segment is the distance between the 

distal oross section of the injection area and a cross section of the substantially 

fibced arteiy segment located at equal distances firom the proximal cross section 

and firom the distal cross section of the substantially fixed arteiy segment The 
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method can furfher comprise the step of creating the model of the substantially 
fixed segment of the artery of the subject The method can further comprise the 
step of determining the projection angle £rom the model, or the step of 
compensating for the non-perpendicularity of the substantially fixed segment of 
the artery of the subject TThe method can furttiOT comprise the step of registering 
the angiograms with the modeL Within the method, the model can be a three- 
dimensional model. The method can further comprise a step of determining TDVH 
grades firom local gray level curves in multiple points of the artery. 

Yet another aspect of the present invention relates to a method for determining 
the arterial reserve of a subject, the method comprising the following steps: 
receiving one or more first models, representing two or more substantially fixed 
segments of one or more arteries of the subject at a plurality of points in time 
associated with and during one or more parts of one or more first heart beat 
cycles; injecting contrast agent into the artery of the subject, said contrast agent is 
delivered to two or more substantially fixed segments being in a non-hyperemic 
state; determining firom one or more first angiograms taken at a predetermined 
projection angle, a first set of density curves, one density curve for each of the 
two or more substantially fixed segment; injecting the subject with substance that 
simulates hyperemia; receiving one or more second models, representing the two 
or more substantially fixed segments of the artery being in a hyperraalc state, at 
times corresponding to a plurality of points in time associated with one or more 
parts of one or more second heart beat cycles; injecting the contrast agent to the 
artery of the subject, said contrast agent is delivered to the two or more 
substantially fixed segments being in a hyperemic state; determining firom one or 
more second angiograms taken at a predetermined projection angle, a second set 
of density curves, one density curve for each of the two or more substantially 
fixed segment The method can further comprise the steps of: determining a first 
time shift between a first density curve taken firom the first set of density curves 
and a second density curve taken firom the second set of density curves, the first 

and the second density curves corresponding to a first substantially fixed segment; 
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detennining a second time shift between a third density curve taken from the first 
set of density curves and a fourth density curve taken from the second set of 
density curves, the third and the fourth density curves corresponding to the second 
substantially fixed segment; and determining the arterial reserve as the ratio 
between the first time shift and the second time shift. The method can fiirther 
comprise the steps of: shrinking the first set of density curves in a shrinking factor 
so that maximal similarity occurs between the shrunk first set and the second set; 
and determining flie arterial reserve as the shrinking factor. Alternatively, the 
method can fiirtho: comprise the steps of: stretching the second set of density 
curves in a stretching factor so that maximal similarity occurs between the first set 
and flie stretched second set; and determining (he arterial reserve as the stretching 
factor. 

Yet another aspect of the disclosed invention relates to an apparatus for 

determining the arterial reserve of a subject from two or more images, the 

apparatus comprises: a corcqponent for receiving a model of one or more 

substantially fixed segment of one or more arteries, the segment having a volume; 

a gray level extraction component for extracting the gray level representing the 

material filling rate and diminishing rate along the at least one arteiy; a deiisity 

curve construction component for constructing a density curve associated with the 

substantially fixed segment from the gray levels; a curve fitting component for 

determining a parameter associated with the density curve or a part thereof; and 

an enhanced artery reserve compon^t for determining the arterial reserve as a 

ratio between a first parameter associated vsdth a first density curve associated 

with the substantially fixed segment being in a hyperemic state and a second 

parameter associated with a second density curve associated the substantially 

fixed segment bemg in a npn-hyperemic state. The apparatus can fiirther comprise 

a segmrat volume component for determining the projection angle or the volume 

of the substantially fixed segment of an at least one artery from the model. The 

apparatus can fiirther comprise a relative arterial reserve component for 

detennining the relative arterial reserve between the first artery and a second 
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artery, said relative arterial rraerve being the ratio between the arterial reserve of 
the first art^ and the arterial reserve of the second artery* The apparatus can 
further comprise one or more image acquiring devices. The apparatus can further 
comprise a device for transferring images acquired by an image acquiring device 
to a processing imit, the processing unit comprises one or more input and one or 
more output devices for receiving input and presenting output to a user. The 
apparatus can further comprise a storage device for storing the images or the 
determined arterial reserve values. 

Yet another aspect of the disclosed invention relates an apparatus for determining 
the arterial flow of a subject from two or more images, the apparatus comprising: 
a component for receiving a model and volumes of one or more substantially 
fixed segment of one or more arteries; a gray level extraction component for 
extracting gray level representing the material filling rate and di m inis h ing rate 
along the artery; a density curve construction component for constructing a 
density curve firom the gray levels; a curve fitting component for fitting a curve to 
the density curve or a part thereof; a fractional contrast material volume 
component for determining the artery segment output during the one or more parts 
of the one or more heart beat cycles of the artery. Tlie apparatus can further 
comprise a segment volume component for determining the volume of one or 
more substantially fixed segments of one or more arteries. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
The present invention will be understood and appreciated more fiilly 
from the following detailed description taken in conjunction with the drawings in 
which: 

Fig. 1 is a schematic illustration of an exemplary environment in which 
one embodiment of the proposed invention is used; 

Fig. 2 is a schematic block diagram of the software or firmware 
computational components of the proposed invmtion, in accordance with a 
preferred embodiment of the disclosed invention; 

Fig. 3 is a schematic flow chart of the proposed mediod for measuring 
the coronary reserve and the relative coronary reserve, in accordance with a 
preferred embodiment of the disclosed invention; 

Fig. 4 is a graph showing the gray level representing the contrast agent 
concentration in a specific location in an artery, in accordance with a preferred 
embodiment of the disclosed invention; ¥ig. 5 is a graph showing the gray level as 
representing the concentration of contrast agent as a function of time and location, 
in accordance with a preferred embodiment of the disclosed invention; 

Fig. 6 is a graph illiistrating the analysis of the graph shown in Fig. S. , 
in accordance with a preferred embodiment of the disclosed invention; 

Fig. 7 is a schematic block diagram of the software or firmware 
computational components of an alternative preferred embodiment of the 
disclosed invention; 

Fig. 8 is a schematic flow chart of the proposed method for measuring 
the coronary reserve and the relative coronary reserve, in accordance with an 
alternative preferred embodiment of the disclosed invention; 

Fig. 9A is an illustration of an artery segment during contrast material 
injection, in accordance with a preferred embodiment of the disclosed invention; 

Fig. 9B is an illustration of an artery segment after a contrast material 

injection, in accordance with a preferred embodiment of the disclosed invention; 
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Fig, 10 is a density, curve, in accordance witli a preferred embodiment 
of the disclosed inv^tion; 

Fig. 1 1 is a density curve of a pulsatile flow, in accordance with a 
preferred embodiment of the disclosed invention; 

Fig. 12A shows an artery with two regions of interest, in accordance 
with a preferred embodiment of the disclosed invention; 

Fig. 12B shows two density curves at a non-hyperemic state, in 
accordance with a preferred embodiment of the disclosed invention; and 

Fig. 12G shows two density curves at a hyperemic state, in accordance 
with a prefened embodim^t of the disclosed invention. 
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DETAILED DESCRIPTION OF IHE PREFERRED EMBODI^^ 

A new and novel apparatus and method for the determination of 

coronary reserve and relative . coronary reserve and o&er flow related 

measurements of a specific coronary artery is disclosed. The proposed metiiod and 

apparatus first determines the 3D model, including the volume of a fixed segment 

of an artery at multiple points of time throughout a heart beat cycle. Alternatively, 

the method uses extrapolation of the 3D model, including the associated volumes 

firom one time segment to another. In yet another alternative, the 3D model and is 

volumes are received firom another source, such as CardiOp by Paieon of Rosh 

Ha'ayin, IsraeL Then, using a contrast agent injected to the subject, the apparatus 

measures the rate at which the injected material fills the segment of the artery 

examined- The filling rate and the volume of the segment of the artery examined 

at the same point in time relatively to the heart beat cycle are used to detennine 

the velocity of the blood flow at that point in time. The velocity of the blood flow 

over the heart beat cycle is measured in two preferred embodiments. In one 

preferred embodiment of the present invention, the velocity is calculated by 

analyzing the gray level through the whole artery segment between two 

consequent images, and then integrating and averaging over one or more fall heart 

beat cycles. In the second preferred embodiment, the velocity is measured using 

multiple points along the artery segment, analyzing the local change over time of 

the gray level- Using the velocity, the process yields the coronaxy flow output at 

resting conditions. Then, the subject is injected with coronary vasodilator sudi as 

adenosine, simulating hyperemia, and the process is generally repeated The 3D 

model and volume determination step can be skipped, since the volume of the 

large arteries mi^t change insignificantly as the result of the injection. 

Additionally, since the heart beat cycle might be shorted after the injection, it is 

possible to down sample the series of volumes. The coronary reserve is then 

calculated as the ratio between the coronary output in hyperemia, and the arterial 

output under resting conditions. This method can be applied to any artery in the 

body of the subject, including the coronary arteries. In accordance with another 
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embodiment of the present invention, in order to determine the relative coronary 

reserve, the whole process is repeated for a second coronaiy artery, which is 

known to be healthy, and a coronary reserve is determined for the second artejy. 

The ratio betwe^ the two coronaiy reserves is the relative coronary reserve. 

Reforing now to Fig. 1 depicting an exemplary environment in which 

one preferred embodiment of the present invention is used. As shown in the figure 

in question, a patient subject 4 is Ijdng down in an operating room, being 

calfaeterized. The catheter 8 is inserted into the subjects body at the groin, or at 

another location and is moved through arteries to Hie area of the heart Any known 

methods of catheterized can be used in conjunction with the present invention. At 

the physician's discretion, the catheter emits a certain amount of contrast agent 

into an artery of the patient. A computerized angiography device 12 directed at the 

tip of the catheter is then activated to produce one or more angiograms of the field 

of view. The angiograms depict the areas vdiere the contrast agent is present in 

various levels of gray color, in accordance with the concentration of the agent 

Tlie resulting images depict the artmes that contain contrast agent Hie higher the 

imaging rate, the better the accuracy of the proposed method. In accordance with 

a preferred embodiment, the angiograms are transferred from the catheterization 

woric station 15, which is a standard component of a catheterization laboratory 

(cath-lab) to a work station 16 which uses the images and data transferred from 

the catheterization work station IS. The work station 16 displays information to 

the operator of the apparatus and sets the parameters for injecting the patient with 

contrast agent during the flow analysis. Once the angiograms have been taken, the 

patient himself does not have to be present at the site. The flow reserve 

calculation can be carried on immediately after the angiograms have been taken or 

at a later time, such as a predetermined time set by the user of the apparatus of the 

present invention. In other preferred embpdim^ts of die present invention the 

catheterization work station IS also comprises the apparatus of the present 

invention as an integrated capability. The work station 16 is preferably a 

computing platform, such as a personal computer, a mainframe computer, or any 
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other type of computing platform that is provisioned with a memory device (not 
shown), a CPU or microprocessor device (not shown), and several I/O ports (not 
shown). Alternatively, the work station 16 can be a DSP chip (not shown), an 
ASIC device (not shown) storing the commands and data necessary to execute the 
methods of the present invention, or the like. The work station 16 can further 
include a storage device (not shown), storing the coronary reserve determination 
application. The coronary reserve determination application is a set of logically 
inter-related computer programs or components of compute components and 
associated data structures that interact to determine the coronary reserve and the 
relative coronary reserve from the angiograms. The angiograms are preferably 
transferred to the work station 16 via a transferring device such as a pre-defined 
I/O port (not shown), DICOM-implementing interface, or analog lines. The 
angiograms are processed by an apparatus of the present invention. The work 
station 16 has input and output devices, preferably a keyboard, a mouse and a 
display where the physician or another stuff member can view or manipulate the 
angiograms and the products of the {plication. Alternatively, work station 16 
delivers the output to another system. 

In an alternative embodiment of the present invention, the coronary 
reserve is calculated offline, based on angiograms that were taken at an earlier 
time and possibly at a different location. For example, where the apparatus of the 
present invention receives angiograms taken from a patient located in a remote 
area, such as in a different building, or a different, city, state, country and then 
transmitted at the same time or later, such as even after the patient has recovered 
or his condition has worsen, via communication lines to the location where the 
apparatus of the present invention is located, then the angiograms can be analyzed 
and processed by die apparatus at such time and in other location, effectively 
allowing the use of the apparatus at any time or location irrespective of the 
location of the patient, or the time at which the angiograms were taken. Thus, no 
interaction with the patient is required beyond taking the angiograms. 
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Refeiring now to Fig. 2, showing the main software or jSimware 
components of the coronary resoye determination application, generally referred 
to as 46. The segment volume and preferred angle component SO is a compute 
program or part of a program, such as CardiOp by Paieon of Rosh Ha'ayin, Israel, 
when used frame by frame or when used for 3D modeling and volume 
extrapolation. The segment volume and preferred angle component 50 takes as 
input a series of images of a fixed segmoit of an artery, taken at 7^ equally spaced 
points in time tfarou^out a heart beat cycle, computes a three-dimensional model 
ofthearteiysegmeat, and from that model determines a series of numbers: 

TO . 

where i denotes the index of the relevant point in time, and is the 

volume of the artery segment at time /, The segmrat volume and preferred angle 
component 50 also outputs an imaging direction in which the artery segment 
should be imaged. This direction minimizes the foreshortening of the artery 
segment of interest, and is therefore substantially perpendicular to the artery 
segment Since the artery segment is generally not strai^t, the three-dimensional 
model of the artery segment uses also for compensating for the local non- 
perpendicularity between the imaging view and the artery segment Alternatively, 

step 50 can be skipped when the three-dimensional model, outputs {^/} and the 

projection angle are provided by an extanal source. 

The fractional blood volume component 54 receives as input a scnes of 

images taken immediately after injecting a contrast agent into the artery. The 

images are taken substantially perpendicularly to the artery, using the angle 

calculated by the segment volume and preferred angle for reserve analysis images 

component 50, so that the average gray level of the internal part of the artery in 

the image represents the amount of contrast agent present at the artery segment at 

that instance. The output of the fractional blood volume calculating component is 

the blood volume that passes through the artery segment between time / -1 and 

time /, for ail / between 1 and n. Hie registration of images taken at the 
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different stages of ilie process is «iabled by using the three-dimensional model X)f 
the artery segment, raeated or received in step 50. Since the 3D model of the 
arteiy and its imaging geometry are available, tiie registration task should be 
performed on lateral shift only, tiiuis minimizing registration errors. 

The laactional velocity component 56, determines the average velocity of the 
blood flow at each time segment i.Two preferred methods for velocity measure 
are detailed hereinaft^. 

In one preferred method associated with the present invaition, the fractional 
velocity component 56 takes as input flie change in the average gray level of the 
pixels depicting the inside of the artery in images taken at the /-I and the i-th 
points of time, representing the velocity of the contrast agent flowing tiirough the 
artery segment at that time (the higher the change in the gray level - tiie higher the 
velocity of the material, and vice versa) . Therefore, the time derivative of the gray 
level represents the velocity of tiie material through tiie artery at tiiat time. Using 
the volume of tiie arteiy segment, Vg, at the same point of time relatively to the 
heart beat cycle, and the velocity of tiie contrast material yields the volume of the 
blood passing in the artery segment during each time slot When usmg this 
embodiment, the amount and rate of tiie injected contrast mataial should not 
cause tiie gray level image to reach saturation. The gray level over all the pixels of 
tiie cross section of tiie artery at location / along tiie artery segment in image i is 
denoted by g;(/) . The projected gray level at every point / along tiie artery 
segment m image i, which takes mto account tiie angle between tiie lihe of sight 
of tiie angiogram and tiie local direction of tiie artery, T , is denoted by 
g (G ,*(/), T ) . For example, the function g can take the form of: 

g(G;(0.^) = Gr(0*sin(T) 
Then, the projected gray leyel over tiie whole artery segment is calculated by the 
formula: 
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L 

Hie difference betweoi the gray levels of two consecutive images, is 
therefore: 

The liaear approximatioa of the time dCTivative of the gray level, 

Msi * 
/Si 

where Af is the duration of the i-th time segment, represents the velocity of the 
contrast agent along the arteiy segment However, to switch from gray level 
dbange to velocity, this ratio should be calibrated using a predetermined 

function J \ > . When tiie contrast agent leaves the artery segment, the 

change in tibie gray level is negative^ However, the velocity of the blood is still 
positive. This sign inversion is also taken care of by the function / . 
Therefore, 

At 

where 5? is the average velocity of the flow in the arteiy segment during the /-th 
time slot . 

Referring now to Fig. 4, showing an alternative embodiment for the 
fractional velocity component 56 of Fig. 2. Fig. 4 depicts the gray level at a fixed 
point X at the arteiy segment after injecting contrast agent, as a function of time. 
On the left hand side 1 10 of the graph of Fig. 4 the gray level depicted is low or 
non existent, representing a state where contrast agent is not present in the 
relevant segment to be examined. When injecting the contrast agrat to the artery, 
the contrast agent flows through the arteries and reaches the relevant segment, 
thus an increase 1 12 in the gray level is depicted in the graph of fig. 4. Next, the 
contrast agent gradually fills the segment until such time when the segment is 
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saturated and a steady state 114 is achieved. At that time, the segment is fully 
colored by the contrast agent. 

Turning now to Fig. 5 showing a family of graphs of the type shown in Fig. 
4, on a common coordinate system. Each graph 116, 117, 118, 119 represents the 
changing gray level relatively to the time, at a different location along the artery 
segment, X. In the example shown, graph 119 depicts the part of the segment of 
the artery to which the contrast agent arrives initially, while graph 1 16 depicts the 
part of the segment to which the contrast agent arrives last. 

Referring now to Fig. 6 showing an illustration of a top view of the three- 
dimensional coUection of graphs 116, 117, 118, 119 of Fig.5, where the gray level 
is presented by a pattern. The higher the gray level at Fig.5, the denser the pattern 
is, at Fig- 6. It should be noted, that although the gray level at Fig. 6 are discrete, 
since Fig. 6 is exemplary only, the gray levels are actually continuous. The time 
immediately after the injection and before the contrast agent reaches a certain 
point in the artery is the np-contrast region 120 in Fig. 6. The transient region 124 
represents the segment location whisre and time in which the contrast agent 
gradually fills the artery. The saturation area 128 represents the segment location 
where and time in which the whole relevant segment is filled with the contrast 
agent The tangent of the angle a 132 between the major axis 136 of the transiCTit 
region 124 and the time axis, represents the velocity of the blood flow. This 
me&od provides the velocity 5/ for the i-th time slot of the heart beat cycle. 
Repeating the injection and the analysis for different time segments will provide 
the velocities { Si }for the full heart beat cycle. 

The time interval between the injection of the contrast agent and the arrival 
of the agent to a certain location along the artery can be used for more accurate 
estimation of TIMI grade, by applying the velocity at multiple points meithod at 
one point, for bolus arrival time measurements. 

Referring now bade to Fig. 2, the average artery segment output calculating 
component receives the Si for the z-th time slot, and calculates the blood volmne 
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flowing through Ae aiteiy segment during the z-th time slot AS/, using the 
formula: 

dBi— 

Multiplying the velocity Si by the duration of the time slot yields the distance 
traveled by the material during the time slot Dividing this ratio by the length of 
the segment, X, yields the relative part of the segment traveled by the nciaterial. 
Multiplymg this quantity by the voiimie of the artery segment at that point in time, 
yields the quantity of blood that flew through the segment during this time slot. 

The average artery segment output component, 58^ first determines the 
overall volume of blood flow through the artery segment during the heart beat 

c/cle,^ . B is determined as the summation of the volumes of the blood flows 

over all time slots: 

The average artery segment output, is the total volume, divided by the duration 

of the heart beat cycle, T. 
Therefore 

H 

If the overall volume of blood flow through the artery segment during 
the heart beat cycle, -8 , or the average artery segment output, Q, are of interest^ 

their values are output by the average artery segment output component, 58. 

The coronary reserve component 62 takes as input the average artery 
segment output in rest conditions, Qr^ and the average artery segment output in 
hyperemia, Qh. The coronary reserve is the ratio: 
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Or 

In a preferred embodiment of the present invention, relative coronary 
reserve calculation component 66 takes as input the CFR of two arteries, CFRa 
which is the artery diseased or suspect of being diseased and CFKb which is of a 
healthy artery and calculates their ratio: 

It will be understood that if a healthy artery is not present than the 
relative coronary flow will not be calculated. 

Referring now to Fig. 7 which shows an alternative to components 54, 
56, 58, and 62 of Fig. 2, and to Fig. 8 which shows an alternative to step 96 of 
Fig. 3, by using the following analysis. Blood does not behave like a Newtonian 
fluid. Blood is a bi--phase material, containing plasma and blood cells, and 
therefore behaves like pseudoplastic material. In addition, due to the Farhaeus- 
Lindqvist effect, blood cells move faster than the plasma in the center of the 
vessel, and since the coronary arteries are not long enough, the velocity profile is 
not fuUy developed. All these factors cause the blood velocity profile to be non- 
compatible with a parabolic model and there is no analytic way to describe all the 
behaviors simultaneously. Therefore, an empiric model for the blood yelocily, 
developed by Verhoeven at 1985 that determines the blood flow is adopted. 
According to the model: 

(1) m = , 

Wherein: 

is the radius of the vessel; 
r is the distance of the measured point firom the center of the vessel; 
Vfn is the maximal velocity, which is the velocity at the center of the 

vessel; and 
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V(r) is the velocity at distance r from the center of the vessel, 

Refening now to Fig. 9A, showing an artery segment during contrast 
material injection, and Fig* 9B showmg the same segment at time t after the 
injection. The contrast material is injected at area 204, denoted by L, and the 
injection point is considered to be the most distal point drea 204, i.e. 221. The 
region of interest (ROI), for which the gray levels are measured resides between 
212 denoted by Xi and 216 denoted by X2. Distance 208 is the mean distance 
between the injection point and the region of interest, and is denoted by X .The 
distance referred to as 200, is equal to 2*iJ. The direction along which Hxc blood 
flows is denoted by x and referred to by arrow 220, which starts at the injection 
pomt 221, On Fig. 9B, the contrast material is shown to be di^ersed as a bolus 
between two borderlines, the upper borderline denoted by 224 and the lower 
borderline denoted by 228. Upper borderline 224 and lower borderline 228 are 
deterniined using the velocity profile shown in (1). 

Using image processing methods, includinjg edge detection and gray 
level analysis performed over angiograms showing a cross section of the artery 
taken at multiple points in time during a heart beat cycle when the contrast 
material is dispersing in the artery segment, and the information related to the 
volume of the artery segment as obtained from the 3D model constructed earlier, 
the density values of the contrast material within the relevant artery segment are 
evaluated. Then a density curve, describing tiie percentage of the contrast material 
out of the total volimie of the segment, is constructed from the density values. The 
density curve is built by density curve construction component 140 of Fig. 7, 
using the gray level analysis of the captured angiograms. 

Geometrical analysis provides the voltraie bounded between upper 
borderline 224 and lower borderline 228 of Fig. 9B. By performing two 
reductions detailed below, the following formula for the volume of the contrast 
material present at time t within the ROI is obtained: 

(2) C^,,jr2(0 = ^Hx2-xl).[l-e'e * ] 
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Wherein the applied reductioiis aie as follows: 

• Dealing only with upper borderline 224 ignoring lower borderline 228. 

• Relating to instead of integrating between A!/ and JTz 

Equation (2) provides an approximated mathematical description of the density 
curve built from the angiograms. Hierefore, Ihe value of the plateau in tiie graph 
is equal to llie constant «R * (x2 -xl) . 

Referring now to Fig- 10, showing a density curve. By using curve fitting to 
the exponential area marked as 300, the e^onential part of equation (2) , i.e., - 

^is obtained. The curve fitting is performed by curve fitting component 144 of 

X 

Fig. 7. The same analysis applies also when the pulsatile nature of the blood flow 
as shown in Fig. 11. Dq)ending on the nature of the curve, it may be beneficial, 
such as when considering the asceading and plateau parts of the d^ity curve 
shown in Fig. 1 1 to perform the curve fitting to a logariliim of the diensity curve, 
instead of to the curve itself. The determination of said ejqwnential part, being 

^ is performed at step 160 of Fig. 8. Once %is obtained, Vr can be obtained 

X X 

for each r, and when integrated over the axws section and avearaged over the 
region of interest, provides the flow through the. artery at instant t, and when 
averaged oVer a full heart beat provides the average artery segment output The 
determination of tiie ratio is performed at step 164 of Fig. 8, by fiactional contrast 
material volume component 148 of Fig. 7. If the flow itself is of no importance, 
step 164 of Fig. 8 can be skipped and fractional contrast material volume 
component 148 of Fig. 7 can be omitted. In this case, when the only required 
mformation is the CFR, the CFR can be directiy determined by 

ywC^MFwiwto) ^ CrossSecHon{hyperenud) 

X : _ . — , which is performed by enhanced coronary 

^^^!!^^*CrossSection(rest) 

X 

reserve component 152 of Fig. 7. If the cross section of the artery does not 
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. • . Vm 

change betwera rest and hyperemia conditions, the ratio between -^z^st 

X 

hyperemia and at rest provide the CFIL The CFR calculation is performed at step 
168 of Fig. 8, which is different from step 104 of Fig. 3, since it determines tiie 
ratio between tibie two rations, rather than the ratio between two volumes, as in 
Fig. 3. As in Fig. 3 above, the relative CFR is determined at step 108 of Fig. 8, as 
the ratio between the CFR of a diseased artery and the CFR of a healthy artery. 

The proposed steps have a number of advantages: The CFR can be found 
without determining!^, by just keeping the same location of ROI and injection 
area; If the flow itself is required, then X can be approximated by the center of 
the ROI; curve fitting is simple as one paramet^ only needs to be evaluated; and 
curve fitting eliminates noise arid provides a more robust solution. 

If the first reduction discussed above, of ignoring the lower borderline is to 
be avoided, the density graph comprises three parts: the ascending part, in which 
contrast material is arriving into the ROI, the part being described by the formula: 


(3) C^ux2i0-K^[l^e'e ^ ] 

the middle part in which material is arriving arid leaving the ROI, and is 
described by the formula: 

(4) C^ux2i0-K^er.[e'-^-^e''] 

and the third part in which the material is leaving the ROI, which is 
described by the formula: 

(5) Cj,^x7(0 = K'e'e'*'- 

Under these conditions, three curves are fitted in order to obtain the three 

parameters: K, — :: — , and — — — 
X x + L 
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Hie disclosed method is not limited to the assmned Verhoeven velocity 

Vm 

profile, which is merely an example, and to the -^parameter. Oth» velocity 

profiles can be used, and dififerent parameters can be obtained from the graphs 
and used for determining the flow output or the arterial reserve, Furfheraiore, the 
method is not limited to one parameter, but can raHier use any number of 
parameters and combinations thereof. 

As yet another alternative to steps 54, 56, 58, and 62 of Fig, 2, a model of 
considering two ROIs, a proximal one and a distal one is disclosed Fig. 12A 
shows a proximal ROI denoted by 3 12 and a distal ROI denoted by 3 16. Fig. 12B 
shows a density curve for proximal ROI 312, marked as 320 and a density ciarve 
for distal ROI 316 marked as 324, both density curves match the respective ROIs 
at resting, i.e., non-hyperemic conditions. Fig. 12C shows density curve 321 
corresponding to proximal ROI 3 12 and density curve 325 corresponding to distal 
ROI 316, both at hyperemic conditions. A method known in the art for 
determining the CFR is by determining the ratio between: 1. the time shift 
between density curves 320 and 324, and 2. The time shift betweai density 
curves 321 and 325, is the CFR. The time shifts can be found, for example by 
cross-correlation between curves 320 and 324, and cross-correlation between 
curves 321 and 325. The disclosed method proposes an improved and more 
accurate alternative, being determining the CFR as the ratio between: 1. the time 
shift between density curves 320 and 321, and 2. the time shift between density 
curves 324 and 325. 

In yet another alternative, the curves at Fig. 12B are shrunk, until maximal 
similarity between the shrunk curves of Fig. 12B and the curves of Fig. 12C is 
achieved. The shrinking factor is the CFR. Alternatively, the curves at Fig. 12C 
are stretched, until maximal similarity between the stretched curveis of Fig. 12C 
and the curves of Fig. 12B is achieved. The maximal similarity can be obtained, 
for example by oross correlation. 
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The proposed methods can be fturther improved by employing a few 
enhancements to the injection procedure. A first enhancement comprises 
performing the injection when tiie flow in the artery is minimal, i.e., during the 
systole (this is true only for the coronary artery, unlike the rest of the body), which 
provides two advantages. The first advantage is that the shape of the contrast 
material bolus is more distinct when there is little or no flow around it, i,e., less 
extenml forces are applied to the material. The second advantage is that the 
contrast agent causes minimal intermpt to the measured flow. Another proposed 
enhancement is to inject the contrast material continuously throu^out an integer 
number of heart beat cycles, i.e., a time period which is a product of an integer 
number and the duration of a heart beat cycles of tiie subject. Using this 
enhancement, the upper and the lower borderlines, as shown in Fig. 9B will be 
synchronized, i.e., they will occur at every point along the artery segment at 
exactiy the same phase of the iieart beat cycle. Yet another enhancement is to 
create a compact bolus of contrast material, by injecting the materisd radially, i.e,, 
sideways towards the walls of the artery rather than from the tip of the catheter in 
the distal direction. The radial injection will significantiy straighten the upper and 
lower borderlines, so that the relevant artery segment will assume a substantially 
cylindrical shape. 

It will be appreciated by persons skilled in the art that various other methods 

and nuances of determining the velocity or the output of the artery segment from 

gray levels exist The methods originate from different formulas representing the 

blood flow, either analytical or empirical, different analysis and reductions and 

different mathematical methods. 

Referring now to Fig. 3, showing the flow diagram of the process of 

calculating the coronary reserve of an artery. In step 84 at least one pair of images 

of the artery are taken at multiple times throughout a heart beat cycle. In each pair, 

the two images are taken from different projection angles. From each pair of 

images, a three-dimensional model of the artery segment is generated, and the 

volume of a fixed part of the artery is determined by a relevant software program 
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or part of a program, such as CardiOp made by Paieon from Rosh Ha'Ayin, Israel. 

Hie part of the arteiy of mtaest is defined, for example, by a fixed distance from 

a stenosis in the distal direction and a possibly different distance in the proximal 

direction. Instead of a stenosis, any other fixed point can be used as a reference 

point for the artery segment Preferably, as long as possible, artery segments with 

no isubstantial branching are chosenL However, branches flow can be subtracted 

for accurate measurement and calculations when the artery segments do include 

branches. The fixed length can be anywhere from about a few millimeters to about 

15 centimeters fit>m the fixed point in the proximal and the distal directions. 

During step 84, the projection angle which is most perpendicular to the artery 

segment is calculated form the three dimensional model of the artery segment. 

Step 84 can be accomplished also by extrapolation of a 3P model performed at 

one point in time during the heart beat cycle to the ^tire heart cycle. 

In step 88, contrast agent is injected to the art^y through the use of the 

catheter as shown in Fig. 1 or throu^ the use of other means such as a needle, 

intravenous tube and the like or a specific controlled injection system. In step 92 

multiple images are taken from the projection angle calculated in step 84, 

throughout a heart beat q/cle, immediately after the contrast agent injection. The 

images are taken at points in time corresponding to the times at which images 

were taken at step 84, i.e,, at the same times relatively to fhe heart beat cycle. The 

more images taken during the heart beat cycle, tihie better is the resolution received 

and therefore the accuracy m providing the CFR and rCFR. The number of images 

taken depends on the duration of the heart beat cycle and on the image-acquiring 

rate of the imaging device. Determining the exact sequence of images taken 

tbroughout exactly one fiiU heart beat cycle is automatically achieved by the 

proposed system, using data from the catheterization work station IS of Fig. 1, or 

any other synchronization equipment. In step 96, tiie average artery segment 

output calculation throughout the heart beat cycle is determined using the 

firactional blood volume calculation coraponent 54 of Fig. 2. and the average 

artery segmrait output calculation 58 of Fig. 2. 
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In Step 100, a coronary vasodilator, sudi as adenosine, is injected to the 
subject, and steps 84, 88, 92, and 96 are tepested, so that the average artery 
segment output is calculated for the arteiy in hyperemia. Optionally, step 84 of 
determining the volume of tiie arteiy segment at multiple pomts can be skipped, 
and the volumes detennined prior to the injection of the coronary vasodilator can 
be used also subsequent to tiie injection, taking into account tiie shortened heart 
beat cycle. This is possible since the coronary vasodilator mainly dilates the small 
vessels that infuse blood into the muscles rather than the large vessels whose 
volume does not change significantly. In addition, since the heart beat cycle is 
shorteaied due to the mjection of coronal vasodilator, it is possible to ignore some 
of the volume values that were collected prior to the injection and use only a 
subset Then, in step 104, the coronary reserve is calculated using the coronary 
reserve calculation component 92 of Fig. 2, by determining the ratio between the 
average artery segment flow output followiiig the adenosine injection (i.e. in 
hyperemia) and the average artery segment flow output prior to the adenosine 
injection (in rest condition) . 

In order to determine the relative coronary reserve, steps 84, 88, 92, 96, 
100, and 104 are repeated for a second artery segment, i.e., the CFR is determmed 
for a second artery segm«it Then, in step 108, the relative coronary reserve is 
determined by the ratio between the coronary reserves of the first and the second 
artery segments. 

Persons skilled in tiie art will appreciate that in order to detMiiime the 
CFR and the relative CFR with higher degree of precision, the described process 
can be performed over multiple heart beat cycles or parts thereof rather than one, 
thus in(xeasing the averagmg accuracy. Since the process averages tiie volume of 
blood flowmg through an arteiy segment oyar time, tiie number of heart beat 
cycles considered in resting conditions and in hyperemia need not be equal. 
Similarly, the number of heart beat cycles considered for the diseased artery and 
for a healthy arteiy need not be equal as welL 
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Hie above shown exaiiq)les serve merely to provide a clear 
und^tanding of tiie invention and not to limit the scope of the present invention 
or the claims appended thraeto. Persons skilled in the art will ^preciate that other 
variants of the metiiod and systems can be used in association with the present 
invention so as to meet the invention's goals. Different methods of determining 
the volume of an artery segment, or of determining the average flow through an 
artery segment can be employed. 

The presented m^od and apparatus are innovative in teams of using 
the volume and the flow information of the artery segment either for a specific 
time slice or throughout a fiill heart beat cycle with timing adjustmrats of the 3D 
model and volumes to the velocity measurements. The proposed invention yields 
the CFR and the relative CFR with higher degree of precision, without requking a 
hi^er degree of invasiveness than a standard catheterization. The invention 
carries out the calculations based solely on an^ograms, and does not require 
additional equipment or special e}q>ertise on the side of the physician, thus it is 
easy and cheap to employ. 

The method is accurate, since the data concerning the volume of the 

artery and the velocity of the blood through the artery are Collected independently, 

thus avoiding the interrelations between the &ctors. Furthermore, since the 

structure of flie artery is found first, it enables the deteaminatiOTi of optimal 

projection view to be used for the velocity determination stage in ord^ to 

minimize the need for unaging conditions compensation. However, iif such 

compensation is required, it is best determined onee the artery's structure and 

orientation is Icnown. The data for the stages is collected form tfie artery only, and 

not form other areas captured in tiie an^ograms, thus minimizing undeshed 

effects. The method is hi^y accurate also since the data is collected and 

analyzed separately for each firame throughout the heart beat cycle, but the total 

results take into account the information collected tiux>u^out tiie cycle. In 

addition, performing the gray level analysis over the whole artery segment, 

minimizes problems of measurement fluctuations. 
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It will be q)preciated by persons iskilled in the art that tiie present 
invention is not limited to what has been particularly shown and described 
hereinabove. Rath» the scope of the present invention is defined only by tiie 
claims which follow. 
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CLAIMS 

What is claimed is: 

L A method for d^ennining the arterial reserve of a subject having a blood 
5 flow, said blood flow having a velocity below or equal to at maximal velocity 
value, the method comprising tiiie following steps: 

receiving an at least one first model, said model representing an at least one 
substantially fixed segment of an at least one artery of ihe subject at a 
plurality of points in time associated with and during an at least one part of an 
10 at least one first heart beat qrcle, said artery s^gmmt having a proxin^ 
section and a distal cross section; 

injecting contrast ag&A into the at least one artery of the subject, said 
subject being in a non-hypetemic state, said contrast agent is injected at an 
injection area having a proximal cross section and having a distance fix>m the 
15 at least one substantially fixed segmmt; 

detennining an at least one first parameter fix>m an at least one first 
angiogram representing the substantially fixed artery segmmt, said first 
angiogram taken fix)m a projection angle; 
iryecting the subject with substance that simulates hyperemia; 
20 receiving an at least one second model, representing the at least one 

substantially fixed segment of the at least one artery of <he subject at a 
plurality of points in time associated with and during an at least one part of an 
at least one second heart beat cycle; 

injecting the contiast agent to the at least one artery of the subject, said 
subject being in a hyp«emic state, said contrast agent is injected at the 
injection area; 

detennining an at least one second parameter from an at least one second 
angiogram representing the substantially fixed segment, said second 
angiogram taken form a projection angle; and 
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determining the artoial resarve as the ratio between one of any of the at 
least one first parameter or a oombination th^eof and any of the at least one 
second parameter or a combination thereof 

2. The method of claim 1 wherein the first parameto: is a ratio between the 
5 maximal velocity value of blood within the artery segment when the subject is 

at a hypcTCTiic state and the distance between the substantially fixed segment 
and the injection area. 

3. The method of claim 1 wherein the second parameter is a ratio between the 
■ maximal* velocity value of blood within the artery segment when the subject is 

10 at a non-hyperemic state and the distance of the artery segment fiom the 
injection area. 

4. The method of claim 1 wherein determination of the at least one first 
parameter ooihprises tiie steps of: 

determining fixjm the at least one first angiogram takCTi at a predrtermin«l 
15 projection angle and the at least one first model, a deiisity curve for the at 

least one artery segment; 

obtaining a velocity profile of the blood flow within the at least one 
artery; 

performing curve fitting for the at least one density curve to determine the 
20 first parameter. 

5. The method of claim 1 wharein determination of the at least one second 
parameter comprises the steps of: . . 

determming &om the at least one second angiogram taken at a 
predetermined projection angle and the at least one second model, a density 
25 curve for the at least one artery segment; 

obtaining a velocity profile of the blood flow within the at least one 
artery; 

performing curve fitting for the at least one density curve to determine the 
second parameta. 
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6. Hie method of claim 1 wherein &e distance between the injection area and 
the at least one substantially fixed s^gmrat is the distance betweai the distal 
<^ss section of the injection area and a cross section of the substantially fixed 
artery segment located at equal distances from the proximal cross section and 

5 fiom the distal cross section of the substantially fixed artery segment 

7. The method of claim 1 further comprising the step of oreating the first or the 
second models of the at least one art^y of the subject 

8. The method of claim 1 wherein the first or the second models are three- 
dimensional models. 

10 9, The method of claim 1 fiirther con5)rising the step of detamining the 
projection angle and the volianes of the at least one substantially fixed 
segment of the at least one artery of tiie subject at a plurality of points in time 
associated with and during the at least one part of the at least one first heart 
beat cycle, fit>m tiie first model 

15 10* The method of claim 1 furttier comprising the step of determining the 
projection angle and the volumes of the at least one substantially fixed 
segment of the at least one artery of the subject at a plurality of points in time 
associated with and during the at least one part of the at least one second heart 
beat cycle, fix)m the second modeL 

20 Il.The method of claim 1 further comprising the step of compensating for the 
non-pespendicularity of the at least one substantially fixed segment of the at 
least one art^ of the subject 
12. The method of claim 1 fiirttier comprising the step of registering the at least 
one first angiogram with the first modeL 

25 13. The method of claim 1 fiirther comprising the step of registering the at least 
one second angiogram with the second model. 

14. The method of claim 1 further comprising a step of determining TIMI grades 
ficom local gray level curves in multiple points of tiie artery. 

15. The method of claim 1 fiirtha: comprising the step of determining a relative 

30 arterial reserve as fiie ratio betwe^ the arterial reserve determined for a first 
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artery segmeat and Ihe arterial reserve determined for a second artery 
segment. 

16. The m^od of claim IS wh^^ein the first artery segment is diseased or is 
suspect as bdng diseased and the second artery segment is healthy, 
s 17. The method of claim 1 wherein the arterial reserve is an arterial coronary 
reserve. 

18. The method of claim 1 wherein the contrast agent injection is performed 
during the sj^tole of the subject 

19. The method of claim 1 wherein flie contrast ag^t injection is performed 
10 continuously tiiroughout an integer numb^ of heart beat cycles of the subject 

20. The method of claim 1 wherein the contrast agent is injected radially. 

21. A method for detCTnining the blood flow output of a subject having a blood 
flow» said blood flow having velocity values below or equal to a maximal 

IS velocity, the method conqirisiiig the following steps: 

receiving at least one model, represmting an at least one substantially fixed 
segment of an at least one artery of the subject at a plurality of points in time 
associated with- and during an at least one part of an at least one heart beat 
cycle; 

20 injecting contrast ag^t into the at least one artery of the subject at an 

injection area having a distance from the at least one subs^tantially fixed 
segmeat; 

detennining from an at least one angiogram taken from a projection angle 
and the at least one model, a density curve for the substantially fixed 
25 segm^t; 

obtaining a velocity profile of the blood flow within the at least one 
substantially fixed segment; 

performing curve fitting for the deaisity curve to determine an at least one 
parameter; 
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substituting the at least oue parameter in the velocity profile to determine 
velocity values; and 

integrating tiie velociiy values over a cross section of the substantially 
fixed segment to obtain the arterial segmeat output values of flie blood flow 
within the at least one substantially fixed segment of the at least one artery. 

22. The method of claim 21 wherein the parameter is ratio betwera the ma x imal 
velocity of blood within the substantially fixed segment and the distance of 
the at least one substantially fixed segmrat firom tiie injection point 

23. The method of claim 21 wherem the distaiice between Ihe injection area and 
the at least one substantially fixed segment is the distance between the distal 
cross section of the injection area and a cross section of the at least one 
substantially fixed artery segment located at equal distances firom the 
proximal cross section and fix)m the distal cross section of the at least one 
substantially fixed artery segment 

24. The method of claim 21 further comprising the step of creating the at least 
one model of the at least one substantially fixed segment of flie at least one 
artery of the subject 

25. The method of claim 21 fiirther conqirising the step of determining the 
projection angle fixim the at least one model. 

26. The mefliod of claim 21 fiirther comprising the step of compensating for the 
non-perpendicularity of the at least one substantially fixed segment of the at 
least one artery of the subject 

27. The mefliod of claim 21 fiirfhCT oon^irising Ihe step of regis 

angiograms with the at least one modeL 

28. The method of claun 21 wherein flie at least one model is a three-dimrasional 

model. 

29. The method of claim 21 fiirther conq)rising a step of detocmining TIMI grades 
from local gray level curvies in multiple points of the artery. 
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30, A method for determining the arterial reserve of a subject, the method 
comprising the following steps: 

receiving an at least one first model, representing at least two substantially 
fixed segments of an at least one artery of the subject at a plurality of points in 
time associated with and during an at least one part of an at least one first 
heart beat cycle; 

injecting contrast agent into the at least one artay of the subject, said 
contaist agmt is delivered to the at least two substantially fixed segments 
being in a non-hyperemic state; 

determining firom an at least one first angiogram takm at a predelemained 
projection angle, a first set of density curves, one density curve for each of 
flie at least two substantially fixed segment; 
injecting the subject with substance that simulates hyperenua; 
receiving at least one second model, representing the at least two 
substantially fixed segments of the at least one artery being in a hyperemic 
state, at times corresponding to a pluraHty of points in time associated with an 
at least ODie part of an at least one second heart beat cycle; 

injecting the contrast agent to the at least one artery of the subject, said 
contrast agmt is delivered to the at least two substantially fixed segments 
being in a hyperemic state; 

determining fix)m an at least one second angiogram taken at a 
predetermined projection angle, a second set of density curves, one density 
curve for each of the at least two substantially fixed segment 

3 1 . The method of claim 30 fiirlher comprising the steps of: 

determining a first time shift between a first density curve taken firom the 
first set of density curves and a second density curve taken firom the second set 
of density curves, the first and flie second density curves corresponding to a 
first substantially fixed segment; 

detmnining a second time shift between a third density curve takm firom 
the first set of density curves and a fourdi density curve tsk&i firom the second 
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of density curves, the third and fbe fourth density curves cotresponding to 
tfie second substantial^ fixed segment; and 

determining Ihe arterial reserve as the ratio between the first time shift and 
the second time shift 
5 32. The method of claim 30 fiirttier comprising the steps of: 

shrinking the first set of density curves in a shrinking factor so that 
rn5iyimal similarity occurs between the shrunk first set and the second set; 
and 

^t ^tTtninmQ the arterial resCTve as the shrinking factor. 
10 33. The method of claim 30 fiirttier conqnising ihe steps of: 

stretching ttie second s^ of density curves in a stretching fiictor so that 
m^yimfll similarity occurs betwem the first set and the str^ched second set; 
and 

detemiining the arterial reserve as the stretdiing &ctor. 

15 

34 An apparatus for determining the arterial resoive of a subject fix>m at least two 
images taken at a projection angle, the apparatus comprises: 

a component for receiving a model of an at least one substantially 
fixed segment of an at least one artery, said segm^t having a volume; 
20 a gray level extraction component for extracting gray level 

representing the material filling rate and diminishing rate along the at least 
one artery; 

a dmsity curve construction component for constructing a dmsity 
curve associated with the at least one substantially fixed segment fix>m the 
25 gray levels; 

a curve fitting component for determining a param^ffl associated with 
the density curve or a part thereof, and 

an enhanced artery reserve component for determining the arterial 

reserve as a ratio between a first parameter associated witii a first density 

30 curve associated with the substantially fixed segment being in a hyp^emic 
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State and a second parameter associated with a second density curve 
associated the substantially fixed s^mmt being' in a non-hyparemic state, 

35. The apparatus of claim 34 further comprising a segment volume component 
for deteraiining the projection angle or the volume of the at least one 
substantially fixed segment of an at least one artery fix>m the at least one 
model. 

36. The apparatus of claim 34 fitrfhtf conqnising a relative arterial reserve 
component for d^ermining the relative artmlal reserve between the first artery 
and a sedond artery, said relative arterial reserve being the ratio between die 
arterial reserve of the first artery and the arterial reserve of the second artery. 

37. The apparatus of claim 34 fiirther comprising at least one image acquiring 
device. 

38. The apparatus of claim 34 furtha: comprising a device for transferriog images 
acquired by an image acquiring device to a processing imit, the processing 
unit comprises an at least one input and output devices for receiving input and 
presenting output to a user. 

39. The apparatus of claim 34 furtha comprising a storage device for storing the 
images or the determined arterial reserve values. 

40. An apparatus for detCTiining the artmal flow of a subject from at least two 
images, the apparatus comprising: 

a component for receiving a model and volumes of an at least one 
substantially fixed segment of an at least one artery; 

a gray level extraction conaponent for extracting gray level 
r^esentmg the material filling rate and dimiciishixig rate along the at least 
oneart^, 

a density curve construction component for constructing a density 
curve firom the gray levels; 

a curve fitting conqioneiit for fitting a curve to the density curve or a 
part thereof ; and 
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a fiBctional contrast matmal volume coiiq)oneat for determining the 
artery segment output during the at least one part of the at least one heart 
beat cycle of the at least one artery. 

41. Hie apparatus of claim 40 further comprising a segment volume con[q)onent 
5 for determining the volume of an at least one substantially fixed segment of 

an at lea£ft one artery. 

42. The method as substantially described in the drawings and specification 
above* ... 

10 

43. The appaiBtm as substantially described in the drawings and specification 
above. 
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